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Abstract — In this paper the problem of finding the sparsest (i.e., 
minimum cardinality) critical fc-tuple including one arbitrarily 
specified measurement is considered. The solution to this problem 
can be used to identify weak points in the measurement set, or 
aid the placement of new meters. The critical £-tuple problem 
is a combinatorial generalization of the critical measurement 
calculation problem. Using topological network observability 
results, this paper proposes an efficient and accurate approximate 
solution procedure for the considered problem based on solving a 
minimum-cut (Min-Cut) problem and enumerating all its optimal 
solutions. It is also shown that the sparsest critical fc-tuple prob- 
lem can be formulated as a mixed integer linear programming 
(MILP) problem. This MILP problem can be solved exactly 
using available solvers such as CPLEX and Gurobi. A detailed 
numerical study is presented to evaluate the efficiency and the 
accuracy of the proposed Min-Cut and MILP calculations. 



Notation 

A subset of transmission lines whose removal 

partitions the network into two disjoint parts. 

The set of all transmission lines. 

The graph of the power network. 

The Jacobian of the measurement function. 

A submatrix of H, consisting of the rows and 

columns in the index sets / and J respectively. 

The I th row of mattix H. 

The complement of an index set /. 

The complement of a singleton set {_/}. 

A large scalar constant treated as "infinity" 

in the MILP procedure. 

The number of measurements (rows of H). 

The number of buses (columns of H). 

A subset of buses (i.e., S C f). 

S is used to define a partition of the network. 

Node weights for bus i. 

Edge weights for transmission line {i,j}. 

Edge weights which also account for the bus 

weights connected to {i,j}. wy = Wy + V; + v ;. 

The set of all buses. 

The n x 1 state vector (phase angles). 

Cut capacity of S. The sum of all weights of 

the edges cut by the partition defined by S. 

The sum of all weights of the buses connected 

by edges which are cut by the partition of S. 

Modified cut capacity, where the edge weights 
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I. Introduction 

A. Critical k-tuples 

A modern SCADA/EMS system relies heavily on the state 
estimator, which estimates the power network states (e.g., the 
phase angles of bus voltages) based on measurements such 
as transmission line power flows, bus power injections and 
bus voltages. An important question related to state estimation 
is whether the network is observable or not; whether the 
states can be uniquely determined based on the available 
measurements. This is a central issue of network observability 
analysis (e.g., fl]]-||9l an d the references therein). While the 
measurements are typically placed so that a power network is 
observable, there exist weak points known as critical measure- 
ments. By definition, if a critical measurement is lost (e.g., 
failure of a meter), then the network becomes unobservable 
(i.e., the states can no longer be uniquely determined). The 
notion of critical measurement also plays an important role in 
another vital power network state estimation function, namely, 
bad data detection (e.g., 0, 0, |T0]]-|[I2| and the references 
therein). Specifically, a bad data detection scheme based on 
measurement residual cannot identify whether a meter is faulty 
or not if the corresponding measurement is critical. A gener- 
alization of the concept of critical measurement is a critical 
fc-tuple, where k is any natural number. A critical £-tuple is 
a set of k measurements such that if all measurements in the 
set are lost then the network becomes unobservable. However, 
losing any subset of p < k measurements would not result 
in the loss of observability. A critical 2-tuple is also referred 
to as a critical set, where bad data can be detected but not 
identified (other terminologies include minimally dependent 
set or bad data group, e.g., lfT3l - lfT5l '). Critical ^-tuples of 
larger cardinalities are also of practical interest, as will be 
explained later. How to compute them is the main topic of 
this paper. 

While critical fc-tuple, network observability, and bad data 
detection are closely related, the mathematical tools to analyze 
them are different. While critical ^-tuples provide the mea- 
surements to remove to render the network unobservable, the 
topic of network observability is the opposite. They include 
checking whether a network is observable or not, and in 
case of an unobservable network which parts of it is still 
observable (i.e., finding observable islands). Likewise, critical 
£-tuple computation and bad data detection are separate issues: 
Essentially, the measurement residual based bad data detection 
theory investigates what detection can be achieved beyond 
the limitations imposed by critical measurements or critical 
^-tuples. These techniques are, consequently, not concerned 
with finding the critical ^-tuples. 

Techniques to identify critical 1 -tuples (i.e., critical mea- 
surements) and critical 2-tuples (i.e., critical sets) are widely 
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known (e.g., 0, 113, OJHED)- In 03 and OS *e calcu- 
lation of critical fc-tuples for k > 2 is considered. However, the 
procedure in |fl8l is efficient only for finding critical ^-tuples 
of lower cardinalities (i.e., k < 3), as will be explained and 
numerically demonstrated later in this paper. The computation 
of critical fc-tuples is inherently computationally intensive, 
because finding critical fc-tuples amounts to a combinatorial 
search, as will be discussed in detail in this paper. 

B. Problem Formulation and its Motivation 

The general setup of this paper is the standard state es- 
timation problem over a linearized DC power flow network 
0, 0. The particular problem considered is to find the 
sparsest (i.e., minimum cardinality) critical ^-tuples involving 
at least one arbitrarily specified measurement. By parameter- 
izing the sparsest critical £-tuple problem with the specified 
measurement, it is possible to examine all weak points in the 
network and not just the weakest point at the boundary of the 
network. The precise description of the considered problem 
is as follows. Let m be the number of measurements in the 
power network, and n be the number of states (i.e., the phase 
angles of the bus voltage phasors). It is assumed that m > n. 
Let H G R'" x " be the Jacobian of the state-to-measurement 
function in a linearized model. Denote H(I,J) as the submatrix 
of H formed by including the rows in an index set / and 
the columns in an index set J. Also denote I and / as the 
complements of / and J, respectively. Then according to, for 
instance (Theorem p. 165), the measurements in an index 
set / form a critical A:- tuple if and only if rank(//(7,/)) <n— 1 
for any j. Here j is the index of an arbitrary reference bus, and 
j denotes the set of all indices except j. The sparsest critical 
fc-tuple problem for a specified measurement i can be written 
as: 

minimize card (I) 
i 

subject to rank(Z/ (/,/)) < n — 1 (1) 
iel 

where card(-) denotes either the cardinality of a set or the 
number of nonzero entries of a vector, depending on the input 
argument. Notice that (Q]) does not explicitly impose the con- 
dition that / cannot contain any strictly proper subsets whose 
removal makes H rank deficient. However, this condition is 
always satisfied at optimality. Problem ([TJ requires a combi- 
natorial search of the rows whose removal makes the H matrix 
rank deficient. In general, no efficient algorithm is available 
to exactly solve (Q}. However, specializing (Q3 to the case of 
power system state estimation results in a significant solution 
efficiency gain because of the special structure of H. The 
demonstration of this is the main contributions of the paper. 
Problem (Q]) is motivated from the following applications: 

1 ) Identifying measurements in small cardinality critical k- 
tuples: While not directly solved as an optimization problem, 
(Q]i is addressed in [18] (Definition 1, even though the term 
"critical set" in fT8l has a different meaning than the one here). 
For any given measurement set, lfl8l finds the measurements 
such that optimal objective value of (Q~|) is less than or equal 
to three. This information is used to determine the set of 
additional measurements to be included, so that the network 



becomes more robust to meter failures. Other, but related, 
meter inclusion problems are also considered in lTT9l . [20]. 
This paper, on the other hand, solves (HJ for all i, regardless 
of the corresponding optimal objective value. These include 
the measurements in the critical fc-tuples of cardinalities less 
than or equal to three as in ifTSll . However, the information of 
the sparsest critical fc-tuples of larger cardinalities can be used 
for a measurement inclusion scheme with a more stringent 
robustness requirement. 

2) Planning of measurement sets: Instead of expanding a 
pre-existing measurement set as in |[T8l - ll20ll . it is possible to 
obtain a cost effective yet meter failure robust measurement set 
by removing appropriate measurements from the full set. The 
solution to (HJ with H corresponding to the full measurement 
set can provide insight into which measurements can be re- 
moved. This measurement removal strategy has the advantage 
that it does not assume any pre-existing measurement set 
which can affect the final measurement set. In Section IV-EI 
a numerical example is presented to demonstrate the potential 
of the above measurement removal scheme. 

3) Cyber-security of power networks: Based on a result in 
||2TI (Corollary Q] in Section III-BI ), problem (Q]) is equivalent 
to another cardinality minimization problem: (f2]) to be de- 
scribed in Section lTl-BI Problem (ffjl arises from cyber-security 
analysis of power networks (e.g., ll2TI - lF25l ). In particular, 
Il24ll -I26l analyze the vulnerability of each measurement i 
using (0, where a malicious attacker inflicts "bad data" in a 
critical £-tuple. In this case, sparsest critical ^-tuples of larger 
cardinalities (i.e., > 3) are of interest because the "bad data" is 
intentional instead of occurring by chance. The solution of iff]) 
can be used to identify the weak points in the measurement 
set in the cyber-security setting. 

Finally, note that for complete safeguard against bad data 
or cyber-attack, the set of all critical ^-tuples (in addition to 
the sparsest ones found by solving (Q3) should be computed. 
However, this would require an enumeration which is not com- 
putationally tractable for realistic applications. The calculation 
of the sparsest critical A:-tuples in (Q]) can identify the network 
vulnerabilities, subject to practical computation constraints. 

C. Contributions and Related Work 

This paper presents two methods to solve (Q]). The first 
method is efficient but suboptimal. It utilizes a sufficient 
condition for critical A:-tuples candidates in 0, Ell . This 
condition is topological. In the setting of this paper, the 
condition states that for any set of transmission lines whose cut 
would separate the network into two disjoint parts, removing 
all line and injection measurements associated with these 
transmission lines would make the network unobservable. 
Using this sufficient condition, a restricted version of the 
sparsest critical A:-tuple problem in (Q]) can be stated as follows. 
If the specified measurement i is a line power flow, then 
the corresponding transmission line must be cut. On the 
other hand, if the specified measurement is a power injection 
at a bus, then one of the incident transmission lines must 
be cut. Then the rest of the transmission lines are cut (or 
not cut) in order to minimize the number of measurements 
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removed, while dividing the network into two parts. This 
cut problem, while being a restricted version of (0, is still 
combinatorial. However, if the injection measurements are 
not directly counted towards the optimization objective (to be 
made precise later), then this modified problem becomes a 
classical minimum cut problem (Min-Cut) (e.g., 11271 ). Min- 
Cut admits scalable solution algorithms (e.g., |28|, |29|). The 
solution to the Min-Cut problem can be used as a suboptimal 
solution to (0. In fact, due to [30|-[32|, it is possible to 
efficiently enumerate all optimal solutions to the Min-Cut 
problem and pick the best available suboptimal solutions to 
(0. This is the idea of the first method of this paper. 

Two previous results are related to the first proposed 
method. As mentioned before, lfl"8l addresses (Q}. In lfT8"l . a 
(non-unique) set of measurements is chosen to be the basic 
measurements. Then the critical fc-tuples containing exactly 
one basic measurement can be identified using a matrix 
factorization approach generalizing the one in [3| (Chapter 
4.5.4). To find critical ^-tuples containing more than one basic 
measurements, a recursive application of the matrix factoriza- 
tion approach for finding critical fc-tuples with only one basic 
measurement is required. For larger k, the recursion becomes 
more expensive as there are ("J possible combinations of p 

basic measurements to be included in the critical ^-tuples, for 

n 

different p < k. To solve for all possible i, in total £ ("J 

p=i 

applications of the matrix factorization procedure are required. 
The computation effort is exponential in terms of network size 
(i.e., the number of buses n). In summary, ifTHl is accurate 
but the procedure is efficient only for a sparse measurement 
set (so that critical fc-tuples of high cardinalities will not be 
encountered). The proposed method in this paper, on the other 
hand, is efficient for solving irrespective of the cardinality 
of the critical &-tuple, because (fl]i is approximately solved 
via Min-Cut. However, as it will be numerically demonstrated 
in Section [V] the accuracy of the proposed method suffers 
when the measurement set becomes sparse. In this sense, [18| 
and the proposed method are complementary to each other. 
Another closely related work is EH . which considers a variant 
of 0. In this variant, the sparsest critical fc-tuple also contains 
at least one measurement. However, instead of being user- 
specified, this measurement is chosen by the optimization to 
find the sparsest non-empty critical A:-tuple. Solving for all 
i leads to the solution to the problem in [21 1 but the converse 
is not true. In addition, [21 1 does not pose their problem as 
defined in this paper. Most importantly, the problem in 11211 is 
posed as a submodular function minimization problem l33ll . 
|34|. While theoretically polynomial-time algorithms exist for 
solving this problem (notably the ellipsoid method [35] and 
more recently |36| of which the complexity is 0(m 8 log(m))), 
no practically efficient algorithms for this class of problem 
have been observed. On the other hand, the Min-Cut problem 
encountered by the proposed method can be solved efficiently 
both in theory and in practice. For example, the complexity 
of 11281 is 0(mn + n 2 log(n)). The practical efficiency will be 
demonstrated by the numerical experiment later in this paper. 

The second proposed method, based on mixed integer linear 
programming (MILP), is exact under a mild assumption, but it 



is less time-efficient. The method is based on the equivalence 
between and 0, to be described in Section III-BI This 
means that can be solved by instead solving 0. Previous 
attempts to approximately solve (0 include, for instance, 1221 
describing an attempt to use matching pursuit (e.g., El), 
and (38) about the application of LASSO (39). The MILP 
formulation, based on [26|, does not admit any polynomial 
time solution algorithms in general. However, there exist good 
MILP solvers such as CPLEX (40| or Gurobi gTJ. The 
major novelty of this second contribution of the paper is the 
combination of [ 2 1 1 and [26|. 

D. Organization of the Paper 

The rest of the paper is organized as follows. In Section HIl 
three known theorems from [2| and [21] are reviewed, and 
a corollary is derived. These theorems form the theoretical 
foundation of this paper. Section [TT] also reviews the Min-Cut 
problem, which is an important part of the proposed algorithm. 
Section [Til] describes the first contribution of the paper: a 
Min-Cut based algorithm which makes use of the topological 
characterization of network observability to find the sparsest 
critical fc-tuples. In Section [IV] the second contribution, the 
exact MILP formulation, is derived with some properties 
discussed. In Section [V] some case studies are presented to 
evaluate the performance of the proposed algorithms. Finally 
Section [Vll concludes the paper. 

II. Technical Background 

This section reviews some known results needed for the 
derivation of the contributions of this paper. Theorems adopted 
from known sources are stated without proof. 

A. A Topological Sufficient Condition for Critical k-tuple 
Candidates 

The first statement is adopted from |0 (Theorem 5). It 
provides a sufficient and necessary condition for network 
observability in terms of spanning trees, which are loop-free 
connected subgraphs of the power network retaining all buses 
but subsets of the transmission lines. 

Theorem 1: A power network is observable if and only 
if there exists a spanning tree with an assignment function, 
mapping from the set of the transmission lines in the spanning 
tree to the set of line power flow and injection measurements of 
the original power network. The assignment function satisfies 
the following properties: 

1) Two distinct spanning tree transmission lines map to two 
distinct measurements. 

2) If the line power flow of a spanning tree transmission 
line is measured, then this transmission line maps to its 
own line measurement under the assignment function. 

3) If the line power flow of a spanning tree transmission 
line is not measured, then the injection measurement of 
one of the two terminal buses of this transmission line 
is the value of the assignment function. 

The following theorem, which is the main theoretical basis 
of this paper, is adopted from Theorem 2 of il2~Ti . It states 
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that if an appropriate choice of measurements are removed, 
then it becomes impossible to form any spanning tree with 
an assignment function defined as in Theorem [JJ Hence, the 
statement provides a sufficient condition for finding candidates 
for critical ^-tuples. 

Theorem 2: Let srf be any set of transmission lines whose 
cut would divide the power network into two disjoint parts. 
Then removing all line power flow measurements in &/ and 
all power injection measurements of the buses connected by 
the lines in would render the power network unobservable. 

Remark 1: Any spanning tree of the network contains at 
least one transmission line in si '. However, under the mea- 
surement removal scheme in Theorem 12 it is impossible to 
define any assignment function in Theorem [JJ for this line. 
Hence the network becomes unobservable. 

Remark 2: While Theorem [2] provides the sets of measure- 
ments whose removal would render the network unobservable, 
these sets are not necessarily critical fc-tuples since their 
subsets might also render the network unobservable. 

Remark 3: The original version of Theorem^ as in |2T1 . is 
more general in that it allows the situations in which s/ divides 
the network into more than two disjoint parts. However, the 
method proposed in this paper cannot exploit the additional 
generality. 

B. Sparsest Critical k-tuple Problem as a Cardinality Mini- 
mization Problem 

The following cardinality minimization problem has been 
studied in power network cyber-security (e.g., ||22| . ll25lD : 

minimize card(//0) 

e (2) 
subject to H(i,:)9 = l 

The following theorem, adopted from Theorem 1 in 12TI . 
establishes that the sparsest critical £-tuple problem in (Q~|i is 
equivalent to (f2j). 

Theorem 3: An index set / is a feasible solution to (Q]) if 
and only if there exists a feasible solution 9 in (fJJ such that 
H(j,:)9 — whenever j I. 

Theorem [3] implies the following statement (proved in Ap- 
pendix) establishing the equivalence between ([TJ and (f2j). 

Corollary 1: The optimization problems in (Q~|i and (f2]l are 
equivalent in that 9* is an optimal solution to (f2]) if and only 
if I* — {j I H(j,:)9* ^ 0} is an optimal solution to ([TJ)- 

C. Min-Cut Problem on an Undirected Graph 

Consider an undirected graph ^ = (Y ,£) where Y and S 
denote the set of nodes and the set of edges respectively, and 
let each edge {i, j} G § be weighted with a scalar wij. Let 
5 C Y be any subset of Y . Define the cut capacity function 

5(5) = Y, w ij such that either (a) i G 5 and j ^ 5 

or (b) i £ 5 and j G 5. 

(3) 

For any two distinct nodes s and t , the s — t Min-Cut problem 
seeks to find a partition of Y into Y = 5 U (Y \ 5) such 



that s and t are in different partitions, and the cut capacity 
is minimized: 

minimize 5(5) 

scr (4) 
subject to s G 5 and t S 

For more detail regarding the Min-Cut problem in (0]), see 
for example ll27l . For efficient solution algorithms, see for 
example [28 1, |29|. The Min-Cut problem is a subproblem to 
be solved in the proposed critical A:-tuple calculation algorithm 
to be described in the next section. 

III. Approximate Critical k-tuple calculation via 
Min-Cut Optimization 

A. A Graph-Oriented Optimization Problem Related to (Q 

The sufficient condition in Theorem[2]provides a topological 
characterization of a subset of the solution candidates of the 
sparsest critical £-tuple problem in (Q~|). This characterization 
leads to a graph-oriented optimization problem which is re- 
lated to, but not exactly the same as, (Q]). The development is 
as follows. Denote the power network as ^ = (Y, where 
Y is the set of all buses and S is the set of all transmission 
lines. Then the set srf in Theorem|2l whose cut would partition 
$ , can be characterized by a bus subset S C Y such that 

stf = {{/,;'} I either (a) i e 5 and j or (b) i g S and j G S} 

(5) 

To describe the number of removed measurements associated 
with S (i.e., sf) according to Theorem |2j the following 
definitions are required. Let wy be the number of meters 
on a transmission line {i,j} £ S, and Vj be the number of 
injection flow meters on a bus j G Y. Then associated with 
S, the number of line power flow measurements to remove is 
5(5) defined in ©. In addition, the number of power injection 
measurements to remove can be defined as 

5j n j (5) = £ vj such that either 

j 

(a) ./ G 5 and 3 i G (Y\S) s.t. {/,./} G £ 
or (b) je(Y\S) and 3 i G 5 s.t. {i, j} G S 

Hence, associated with 5, the total number of measurements 
to be removed is 5(5) + 5; n j(5). Lastly, the constraint in (Q]) 
that one specified measurement must be included in the critical 
A:-tuple should be enforced. For simplicity of discussion, for 
the moment it is assumed that the specified measurement is a 
line power flow. The case of power injection will be handled 
in the end of this section. Now suppose the specified line 
power flow meter i is on transmission line {s,t}, then the 
corresponding topological constraint is that s G 5 and t £ S. In 
summary, the graph oriented optimization problem, set up as 
an approximation to ([TJ, is described as: 

minimize 5 (5) + 5i n j (5) 

scr . (6) 

subject to s G 5 and t 5 

Strictly speaking an optimal solution to © is a set of buses, 
with the corresponding set of "cut" transmission lines defined 
in (|5]l. However, it is more convenient to treat the optimal 
solution as the corresponding set of measurements to be 



5 




3 3 



■ meter (] bus 

Fig. 1: Illustration of the modified cut capacity function. Left: 
Metering scenario in (0. Right: Transmission line metering 
scenario pertaining to 5(5) in (0. The number in the paren- 
thesis indicates the number of times a meter is repeated. If 
both 1-2 and 1-3 are cut, the true cost 5 (5) + 5nj (5) is 4. 
However, the approximate cost 5(5) is 5. 

removed, as prescribed by Theorem [2] Solving (0 yields a 
sparse set measurements whose removal makes the network 
unobservable. The numerical experiment in Section [V] will 
demonstrate the usefulness of (0. However, it should be 
emphasized that an optimal solution to (0 is not necessarily 
a sparsest critical £-tuple in 0. The reason is twofold. First, 
since Theorem [2] is a sufficient condition, (0 searches only 
through a subset of the sets of measurements whose removal 
would render the network unobservable. Second, as pointed 
out in Remark |2] an optimal solution to (0 does not even need 
to be a critical A:-tuple. These restrictions will be demonstrated 
by the numerical experiment in Section [V] 

B. Min-Cut Approximate Solution Procedure for (0 

Compared with the tractable Min-Cut problem in (0, (0 
is a combinatorial optimization problem because of the addi- 
tional term 5i n j(5) in the objective function. To overcome the 
computational difficulty, it is proposed in this paper that 5nj (5) 
is indirectly accounted for by solving the following Min-Cut 
problem: 

minimize 5(5) 

Scf (7) 
subject to s 6 5 and t ^ 5, 

where 5(5) is defined according to (0 with modified edge 
weights Wij = Wij + Vj + Vj for all {/, j} G S, with wy, v; and 
Vj defined in Section IIII-AI 5 (5) corresponds to a modified 
metering scenario where an injection meter of a bus is moved 
to all incident transmission line(s). However, this modification 
can lead to overcounting of injection meters as opposed to 
solving (0. See Fig.0for an illustration. Since (0 and (0 have 
the same constraint, the optimal solution to (0 is a suboptimal 
solution to (0. In addition, by the results in [30|-[32|, it is 
possible to efficiently enumerate all optimal solutions to the 
Min-Cut problem in (0. Hence, the best available suboptimal 
solution to (0 can be chosen. However, it is emphasized 
that the strategy of solving (0 can only provide a sparse set 
of measurements including {s,t}, and the removal of these 
measurements makes the network unobservable. From (0 to 
(0 and then to (0 these two transitions induce their respective 
limitations. As it was explained earlier, (0 and (0 are two 



different optimization problems. Moreover, solving (0 is not 
equivalent to solving (0, as illustrated in Fig. The quality 
of approximately solving (0 via (0 depends on the ratio 
between the number of transmission line measurements and 
bus injection measurements. In the extreme case where there is 
no injection measurement, (0 is the same as (0. The accuracy 
and efficiency of the proposed Min-Cut procedure will be 
numerically assessed in Section [V] The following algorithm 
summarizes the Min-Cut based approximate solution proce- 
dure for 0, where the specified measurement is a line power 
flow on a transmission line. 

Algorithm 1: Min-Cut procedure for transmission line 
case: 

Step 1 

In the power network graph & = ('V define arc 
weights Wij as the number of meters on a transmis- 
sion line {i,j} G S, plus the number of meters on 
buses i and j. 

Step 2 

Suppose the specified transmission line is {s,t} 6 S . 
Setup a s - 1 Min-Cut problem as in (0. Solve (0 
using algorithms such as |28l . ||29ll for an optimal 
solution, which is a set of "cut" transmission lines. 
The line power flows measurements and injections 
measurements at the terminal buses constitute a 
suboptimal solution to the sparest critical fc-tuple 
problem in (0. 

Step 3 

Use the results in Il30l - ll32ll to enumerate all optimal 
solutions to (0. Pick the best suboptimal solution to 
among all optimal solutions to (0. 

Even the best available suboptimal solution to might not 
be a critical £-tuple in that there might be a strictly proper 
subset whose removal makes the network unobservable. To 
make sure a critical £-tuple is obtained, an enumeration is 
required to see which measurements in the suboptimal solution 
can be eliminated. However, since the suboptimal solution 
typically contains very few measurements, the enumeration 
is not expensive. 

In the case where the specified measurement i in (0 is 
a power injection on a bus, the following procedure can be 
applied: 

Algorithm 2: Min-Cut procedure for bus injection case: 

Step 1 

Let 5f = ("f ',<?) be the power network graph and 
let i 6 Y be the bus with the considered injection 
measurement. For each j G f such that {i,j} £ S, 
apply Algorithm 1 on transmission line {i,j}. 

Step 2 

Among all solutions provided by Algorithm 1 applied 
to {i, j}, pick the one with the minimum cost in (0 
as the best available solution to 0. 

The numerical examples in Section [V] illustrate the perfor- 
mance of these algorithms. 
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IV. Exact Sparsest Critical k-tuple Problem 
Formulation as a MILP Problem 

Theorem [3] and Corollary Q] state that the sparsest critical 
fc-tuple problem in (Q} can be solved by solving the cardinality 
minimization problem in ([2j. Problem (O can be formulated 
as a MILP problem, as mentioned in ||26l . The key to the 
formulation is the counting of the cardinality of vector H9. To 
achieve this, an additional binary decision vector y £ {0, l} m 
and a scalar constant M > are needed. If M is large enough, 
then the constraint 

\H(j,:)0\<My(j) 

provides a cardinality counting mechanism via y. If \H(j,: 
)8\ >0, theny(j) = 1. If \H(j,:)8\=0, then y(j) can be either 
or 1. However, since <J2J seeks to minimize the cardinality 
of H9, as it will clear shortly, y(j) must be zero at optimality. 
The constant M must be chosen large enough so that it is 
larger than ma\{\H(j,:)9*\} for at least one optimal solution 

9* of (|2). In the special case where all line power flows are 
measured, the method in ll26l can be used to compute M. In 
other cases, the general guideline is that M should be as large 
as possible, before the optimization solver complains about 
numerical difficulties. Suppose 0° is a typical state vector 
under normal operation, then a ^.^g,, ! with some a > 1 
can be a reasonable guess for M. The choice of M is the 
only heuristic part of the otherwise exact sparsest critical k- 
tuple problem formulation. The MILP formulation of (O is as 
follows: 

minimize T.y(j) 

9,y j 

subject to H9 < My 

H9 < My ( g ) 

H(i,:)9 = 1 

y(j) E {0,1} Vj 

Note that since the objective function is Ly(j), at optimality 

for any j such that \H(j,:)9\ — 0, the corresponding y(j) 
must be zero. Hence, =caid(H9). Finally, notice that 

if the measurements in a certain set are considered very 
reliable and are immune from faults, then (0 can be modified 
accordingly by adding the constraint y(j) = for all j E 3P. 

V. Case Study 

Numerical experiment results are demonstrated in this sec- 
tion. All computations are performed on a laptop with an 
Intel Core i5 2.53GHz CPU and 4GB of memory. All Min- 
Cut problems are solved in MATLAB using l42l . which calls 
the libraries from [ 4-3 1 . All MILP problems are solved in 
MATLAB using Gurobi ED via Gurobi Mex 04). 

A. Comparison with the Procedure in [18] 

First, problem (Q~|i, for each possible specified measurement 
i, is solved using three methods. The first method is the 
recursive critical fc-tuple calculation procedure in [18| im- 
plemented by the authors. As in lfT8l . only critical ^-tuples 
containing three basic measurements are sought. The second 



method is the proposed Min-Cut procedure in Algorithm 1 
and Algorithm 2 in Section [III] The third method is the 
MILP procedure (with M = 100) in Section HVl The solution 
of the MILP procedure is used as a reference for accuracy. 
The IEEE 14-bus benchmark system is analyzed, with two 
different measurement sets. The first measurement set is from 
ifTSl (Section IV, measurement set of Scenario 1), containing 
measurements from 9 out of 14 buses and 6 out of 19 trans- 
mission lines. The second measurement set contains all bus 
and transmission line measurements, which may be of interest 
for meter placement. For each i, the procedure in lfl8l (with 
three basic measurements) and Min-Cut might only provide 
an overestimate of the cardinality of the sparsest critical k- 
tuple. Table U lists the percentages of i with overestimation, 
the average overestimation (over all z) and the average relative 
overestimation (relative to the cardinality of the corresponding 
sparsest critical £-tuple). 

Table Q] confirms the statement in Section II-CI that the 
procedure in lfl8l and Min-Cut should be used in different 
measurement settings. [18] performs better in a sparse mea- 
surement set, while Min-Cut is a better choice in a dense 
measurement set. 



TABLE I: Comparison between the procedure in 1181 . Min- 
Cut and MILP for the IEEE 14-bus system 



measurement set 1 (sparse metering) 


|18'| 


Min-Cut 


MILP 


solve time (s) 


0.04 


0.02 


3.6 


percent of meas. with overestimation (%) 





93 





average overestimation 





1.07 





average relative overestimation (%) 





75 





measurement set 2 (full metering) 


["18"] 


Min-Cut 


MILP 


solve time (s) 


10 


0.03 


17 


percent of meas. with overestimation (%) 


26.5 








average overestimation 


4.1 








average relative overestimation (%) 


67.7 









B. The Effect of the Proportion of Line Power Flow Measure- 
ments on the Min-Cut Procedure 

As explained in Section Hill the Min-Cut procedure for (Q~|) 
achieves computation efficiency by approximately counting 
the injection measurements. Hence, the relative ratio between 
the line power flow and bus injection measurements affects 
the approximation quality of the Min-Cut procedure. In this 
subsection, the relationship between approximation quality 
and the proportion of transmission line measurements in the 
network is considered. 

The IEEE 14-bus, 57-bus and 118-bus benchmark systems 
are considered. The network topologies are from MATPOWER 
|45|. For each system, 11 different measurement sets are 
considered. Each measurement set contains all injection mea- 
surements, but the proportions of removed line power flow 
measurements increase as 0%, 10%, . . ., 100%. The removed 
line measurements are randomly chosen. A study similar to 
the one in Section [V-AI is performed, testing only the proposed 
Min-Cut and MILP procedures. The above study can be con- 
sidered as "one sample" of a random experiment involving a 
sequence of 1 1 measurement sets for each benchmark system. 
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The randomness stems from choice of the removed line flow 
measurements. To examine the typical phenomena, the above 
random experiment is repeated five times. Table HI1 shows the 
mean value (over 5 experiments) of the performance and error 
statistics similar to those in Table [] 

While not seen from Table ITT1 the computation time for Min- 
Cut remains roughly the same. The increase in solve time 
ratio (up to about 0.36 when 100% of line measurements 
are removed) is due to the decrease in solve time of the 
MILP procedure. In general, Min-Cut is more efficient than 
MILP. In terms of approximation error, for up to 90% of 
line measurement removal, Min-Cut results in at most 7% of 
measurements whose sparsest critical A:-tuple cardinalities are 
overestimated (the number is down to 3% for up to 40% of 
line measurement removal). On average, the overestimation is 
by about 1 measurement with the maximum observed in the 
experiment being 3 measurements (not shown in Table [TTJ . 
Finally, the average relative overestimation is at worst 35% 
(e.g., overestimation by 1 measurement for a critical 3-tuple). 

C. The Effect of the Proportion of Injection Measurements and 
Arbitrary Measurements on the Min-Cut Procedure 

In this subsection, the experiment in Section lV-Bl is repeated 
for the 118-bus benchmark system with a difference in the 
definition of the measurement sets. Two cases are considered. 
In the first case, each measurement set contains all line power 
flow measurements and different proportions of the injection 
measurements are randomly removed. In the second case, 
different proportions of arbitrary measurements (injection or 
line) are randomly removed in a way that the resulted network 
is still observable. Table [Til] lists the relevant statistics for 
the first case and suggests (from the fourth row) that the 
Min-Cut procedure is more accurate when the transmission 
lines are more densely metered. On the other hand, Table [IV] 
lists the statistics for the second case. In this case, at most 
1 - 118/(118 + 186) w 41% of the measurements can be 
removed before the network become unobservable. However, 
random removal of 40% of the measurements typically results 
in a unobservable network, and hence the corresponding result 
is not shown in Table [IV] Table [IV] confirms again that the 
Min-Cut procedure is efficient and accurate for relatively dense 
measurement sets. 

D. Time Efficiency of Min-Cut and MILP for Large Networks 

This numerical study investigates the possible advantage of 
the proposed Min-Cut procedure for the sparsest critical k- 
tuple analysis for larger scale power networks. The networks 
considered are the IEEE 118-bus, IEEE 300-bus and the 
Polish 2383-bus systems. The topologies of these networks 
are obtained using MATPOWER |45|. On each transmission 
line of the benchmark systems, there are two line power flow 
meters (one from each terminal bus). In addition, all power 
injections are measured. For each of the benchmark system, 
problem $Q is solved using the Min-Cut procedure for all 
possible specified measurement ;. For the 118 and 300 bus 
cases, the Min-Cut procedure is experimentally found to be 
exact, compared with solving the MILP formulation in (|8). 



For the 2383-bus case, only 14 instances of (O are solved. 
For all these 14 instances Min-Cut also provides the correct 
estimates. The computation times for using the proposed Min- 
Cut procedure and solving ([8]l are listed in Table [V] This 
numerical study again assures the efficiency and accuracy 
of the proposed Min-Cut procedure in the case when all 
transmission lines and buses are metered. 

E. Using Critical k-tuple Information for Meter Placement 

While not the main focus of this paper, a possible use of 
the critical £-tuple information in ([TJ is for meter placement. 
Problem ([TJ is solved for all specified measurement i for the 
IEEE 6-bus benchmark system (see Fig. |2). It is assumed 
that the network is fully metered (with 6 injection and 11 
line power flow measurements). The Min-Cut procedure (i.e., 
Algorithm 1 and Algorithm 2) is used to solve (|T). In total, 17 
critical fc-tuples are found. Then for each measurement i, Fig. [3] 
shows the number of critical fc-tuples containing i. However, 
note that by solving (Tj]i for all specified measurement i, all 
sparsest critical fc-tuples are not found. Hence, Fig. [3] only 
shows the lower bounds for the true number of critical ^-tuples 
containing i. Fig. [3] indicates that measurement 12 is probably 
not important (as far as network observability is concerned), 
because only one critical fc-tuple contains it and none of the 
other 16 critical k- tuples will be affected by the removal of 
measurement 12. This is consistent with the network topology 
in Fig. [2] since measurement 12 is the line power flow 
measurement between bus 2 and bus 5 (i.e., the two buses 
with the largest degree). On the other hand, Fig. [3] suggests 
that measurements 2 and 5 (i.e., power injections at bus 2 
and bus 5) are definitely important because they are involved 
in many critical ^-tuples. This is again consistent with the 
topology in Fig. [2] since each of these injection measurements 
can substitute one of the five line measurements in case any 
one of them fails. Finally, note that the same analysis here can 
be carried out for larger scale networks where it becomes less 
obvious from the topology which measurements are important 
or unimportant. 



1 2 3 




4 5 6 

□ bus 



Fig. 2: IEEE 6-bus system. 



VI. Conclusion 

In conclusion, a version of the sparsest critical £-tuple 
problem is considered. The sparsest critical fc-tuple is sought 
for one arbitrarily specified measurement. It is possible to 
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TABLE II: Ensemble mean of the solve time and error statistics for the case study in Section IV-BI 



14-bus 


line meas. removal (relative to total lines) (%) 
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average relative overestimation (%) 
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line meas. removal (relative to total lines) (%) 
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line meas. removal (relative to total meas.) (%) 
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12 


18 


24 


31 


37 


43 


49 


55 
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solve time ratio (Min-Cut/MILP) xl(r 4 


2.1 


2.1 


1.7 


2.0 


2.1 


2.4 


2.7 


3.9 


4.6 


11 


1116 


percent of meas. with overestimation (%) 





0.77 


0.90 


1.4 


2.7 


4.2 


5.4 


6.7 


6.7 


5.7 


87 


average overestimation 
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1.2 
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1.2 


1.2 


1.2 


1.1 


average relative overestimation (%) 





13 


15 


14 


17 


19 


22 


24 
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35 
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TABLE III: Ensemble mean of the solve time and error statistics for the case study in Section IV-CI with varying injection 
measurement proportion 



bus meas. removal (relative to total buses) (%) 
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20 


30 


40 


50 


60 


70 


80 


90 
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bus meas. removal (relative to total meas.) (%) 
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12 


16 


19 


23 


27 


31 


34 


39 


solve time ratio (Min-Cut/MILP) xl0~ 4 


2.1 


3.2 


6.4 


9.7 


18 


33 


61 


108 
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404 


1211 


percent of meas. with overestimation (%) 





0.61 


1.4 


1.4 


1.5 
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2.6 


2.3 


1.4 


1.1 





average overestimation 





1.4 


1.1 


1.2 
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1.2 


1.1 


1.2 


1.1 


0.8 





average relative overestimation (%) 
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14 
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18 


19 


20 


22 
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TABLE IV: Ensemble mean of the solve time and error statistics for the case study in Section IV-CI with varying arbitrary 
measurement proportion 



meas. removal (relative to total meas.) (%) 





10 


20 


30 


40 


50 


solve time ratio (Min-Cut/MILP) xl0~ 4 


2.1 


3.3 


8.9 


21 


110 


345 


percent of meas. with overestimation (%) 





1.0 


5.0 


7.6 


20 


35 


average overestimation 





0.9 


1.0 


1.1 


1.1 


1.4 


average relative overestimation (%) 





12 


24 


28 


45 


69 




8 10 12 

measurement index 



Fig. 3: Number of critical ^-tuples containing any specific 
measurement. 



TABLE V: CPU times for solving all instances of (HJ in the 
118, 300 and 2383 bus systems 



Method 


118-bus 


300-bus 


2383-bus 


MILP 


763 sec 


6708 sec 


(projected) 5.7 days 


Min-Cut 


0.30 sec 


1 sec 


31 sec 



paper demonstrates that the studied sparsest critical fc-tuple 
problem can be formulated as a MILP problem so that power- 
ful MILP solvers such as CPLEX and Gurobi can be utilized. 
On the other hand, by using topological network observability 
results in (0, lETT l. a Min-Cut based approximate solution 
procedure can be derived. The numerical experiment in this 
paper reveals that the Min-Cut procedure is highly accurate 
and efficient when there are a significant number of line power 
flow measurements in the power network. Consequently, Min- 
Cut should be the first method to attempt (over MILP) in this 
scenario. 



identify the weak points in the power network by listing 
all measurements which might form critical fc-tuples with 
small cardinality, even though this is short of a complete 
enumeration of all possible sparsest critical ^-tuples. This 



Appendix 

Proof of Corollary [7] 

Suppose 0* is an optimal solution to (0 and /* is such 
that H(j,:)9* ^ if and only if ; e /*. Then Theorem [3] 
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states that 7* is a feasible solution to _Q) with the objective 
value card(/*) — card(770*). Now suppose 7* is not optimal 
and there exists another feasible solution 7 of (Q_ such that 
card(7) < card(/*). Then Theorem [3] states that there exists 
9, feasible in ©, such that card(77f5) < card(7) < card(/*) = 
card(770*). This contradicts the assumption that 9* is an 
optimal solution to (0. Hence, 7 does not exist and 7* is an 
optimal solution to (Q~|). To establish the converse, suppose 7* 
is an optimal solution to __). Then Theorem [3] states that there 
exists 9* feasible to __> such that card(770*) < card(7*). In 
fact, card(770*) = card (7*) and 9* is optimal to __>. If this is 
not true, Theorem [3] implies that 7* would not be optimal. 
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